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ABSTRACT 

Methods for quantifying environmentally induced decoherence in quantum systems are 
investigated. We formulate criteria for measuring the degree of decoherence and consider 
several representative examples, including a spin interacting with the modes of a bosonic, 
e.g., phonon, bath. We formulate an approach based on the operator norm measuring the 
deviation of the actual density matrix from the ideal one which would describe the system 
without environmental interactions. 
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1 Introduction 



Evolution of a quantum system exposed to an environment is described by the density ma- 
trix and deviates from the ideal, usually pure-state, dynamics. In this work, we consider the 
problem of quantifying this deviation by single numerical measures, derived from the full set 



of the density matrix elements. Depending on their strength and nature, environmental inter- 
actions can lead to various relaxation and even measurement-type processes. We will use the 
term "decoherence" generically. Our focus will be on the environmental effects that represent 
"noise" and cause small deviations from the desired isolated-system dynamics. Establishment 
of the threshold criteria for fault-tolerant quantum computation, 1 has made estimation of the 
noise due to the environmental effects important for evaluating quantum computing systems. 

Superposition of quantum states is crucial in utilizing quantum parallelism for quantum 
computation. Therefore quantum algorithms usually deal with pure or nearly pure states. 
There are actually two possible measures of decoherence: deviation from a pure state, and de- 
viation from a particular pure state, or nearly-pure mixed state, which can be time-dependent. 
Even if the state of the system remains pure, it might deviate from the dynamics desired for a 
particular controlled quantum process. The ultimate goal of the studies of the type reported 
here would be to identify simple numerical measures of degree of decoherence occurring during 
the "clock" times of quantum computing gate functions, to compare with the fault-tolerance 
requirements in various quantum error correction schemes. 

Let us assume that a quantum system of interest, denoted by S, is prepared at time t = in 
some pure quantum state \ip). At times t > 0, the system can be subject to noise, i.e., coupling 
to the modes of the environment, but also to "controlling" interactions required for quantum 
computing. Here we consider a system only exposed to the environmental noise, because the 
degree of decoherence can be usually approximately evaluated for time scales shorter than or 
comparable to the times of quantum control. The quantum system can then be described by 
the reduced density operator, p(t), which can be obtained from the overall density operator 
by tracing over all the environmental degrees of freedom. 

For the initial state, the density matrix has the form p(0) = \<p)(<p\- The interaction 
with the surrounding, which is usually assumed to be a large macroscopic system in thermal 
equilibrium at temperature T, leads to thermalization of the quantum system. The reduced 
density matrix of the system for large times should approach 

p^e^ s /Tr 5 (e- /3Hs ), ast^oo, (1) 

where (3 = 1/kT. Here H s is the Hamiltonian of the system. In the energy eigenbasis, 
the density matrix elements behave as follows, 

Pmn = {<Pm\p\<Pn) ~* 0, as t ► OO (m ^ 7l), (2) 

Pun = (<t>nWn) ~> e"^/ £ e"^* , as t > OO. (3) 

k 

Usually, thermalization sets in as the slowest relaxation process, and its description requires 
a phenomenological Markovian assumption. 2-4 



In this work, we consider several approaches to measuring the degree of decoherence due to 
interactions with the environment. In Section 2, we discuss the approach based on the asymp- 
totic relaxation time scales. The entropy and idempotency-defect measures are addressed in 
Section 3. The fidelity measure of decoherence is discussed in Section 4. Next, in Section 5, we 
present our results on the operator norm measures of decoherence. Section 6 presents a dis- 
cussion of an approach to eliminate the initial-state dependence of the decoherence measures, 
as well as of extensivity properties for multiqubit systems. 



2 Relaxation Time Scales 

Markovian approximation schemes typically yield exponential approach to the limiting 
values of the density matrix elements for large times. 2-4 For the two-state system, this defines 
the thermalization time scales Ti and T 2 , associated, respectively, with the approach by the 
diagonal (thermalization) and off-diagonal (dephasing, decoherence) density-matrix elements 
to their limiting values; see (2,3). More generally, for large times we expect 

Pnn(t) -Pnn(00) OC e~ t/T ™ , (4) 

p nm (t) oc e~ t/Tnm (n^m). (5) 

The shortest time among T nn is often identified as Ti. Similarly, T 2 can be defined as the 
shortest time among T n ^ m . These definitions yield the characteristic times of thermalization 
and decoherence. 

For decoherence and thermalization times, the following inequality commonly holds, T 2 < 
Ti, 3 though this relation usually does not apply within Markovian approximations. Therefore, 
the decoherence time is a more crucial parameter for quantum computing considerations. 
The time scale T 2 is compared to the "clock" times of quantum control, i.e., the quantum gate 
functions, T g , in order to ensure the fault-tolerant error correction criterion T g /T 2 < O (10~ 4 ). 5 

The disadvantages of this type of analysis are that the exponential behavior of the density 
matrix elements in the energy basis is applicable only for large times, whereas for quantum 
computing applications, the short-time behavior is usually relevant. 6 Moreover, while the 
energy basis is natural for large times, the choice of the preferred basis is not obvious for short 
and intermediate times. 6 ' 7 Therefore, the time scales Ti and T 2 have limited applicability in 
evaluating quantum computing error correction criteria. Their obvious advantage is in that 
of all the measures discussed in this article, they are the only experimentally observable time 
scales, as long as no multi-qubit quantum computer was built. 



3 Quantum Entropy 



An alternative approach is to calculate the entropy 8 of the system, 

S(t) = -Tr(plnp), (6) 

or the idempotency defect, called also the first order entropy, 9-11 

s(t) = 1 - Tr (p 2 ) . (7) 

Both expressions are basis independent, have a minimum at pure states and effectively describe 
the degree of the state's "purity." Any deviation from the pure state leads to the deviation 
from the minimal values, 0, for both measures, 

S pure statc(^) Spurc statc(^) 0. (8) 

Let us consider a simple example of a two-level system in the excited state |1) at t — 0, 
which decays to a ground state |0) via the interaction with a reservoir at zero temperature 
T = 0, e.g., the spontaneous decay of an atom in vacuum due to the electromagnetic field. 
Except for very short times, 3 the dynamics of the system is well described by the following 
density matrix, 

Pll (t) = e- rt , p 00 (t) = l-e- rt , (9) 
PoiW = Pio(0 = Poi(0) = pio(0) = 0, (10) 

where T is the rate of the spontaneous decay. 

In Figure 1, we plot the two entropy measures of decoherence as functions of time, for 
T = 1. Both curves illustrate that these measures determine deviations from "purity," but are 
not sensitive to deviations from any particular pure (or mixed) state. It is likely that entropy 
measures are useful in situations when the system has evolved little from a particular initial 
pure state under its own dynamics and environmental influences. Entropy deviations can thus 
be used to quantify degree of decoherence/relaxation for short times, 6 when the system is still 
sufficiently far from any but the ideal pure state (corresponding to evolution under H$ only). 

One can derive 9 the short-time perturbative expansion for s(t). Let R(t) denote the density 
matrix of the system plus environment, and, for definiteness, assume that the system is initially 
in the pure state \ip), unentangled with the environmental bath (B) of modes which are 
described by the density matrix 9. Then, 



p(t) = Tr B R(t), 



(11) 



12 3 

t 



Figure 1: Entropy measures of decoherence for the two-state spontaneous decay model. The 
upper curve: the full entropy. The lower curve: the first-order entropy. 



R(t) = e- iHt R(0)e iHt , (12) 

fl(O) = (K>)M)®0, (13) 
H = H S + H B + V, (14) 

where H is the total Hamiltonian, H$ is the Hamiltonian of the quantum system, Hb is the 
Hamiltonian of the environmental bath, V is the interaction Hamiltonian. Here we use the 
units h — 1, and we omit the direct product symbol in what follows. We get the expansion of 
the form 

s(t) = (t/r s ) 2 + o(f), (15) 



l/(2r s 2 ) = (((V) B ) 2 )s + (((V)s) 2 )b - ((V 2 ) s )b - (((V) s )b) 2 , 



(16) 



where (...)b, (•••)5 denote averaging with respect to 6 1 and p(0), respectively. Here, as usual, 
average of an operator stands for the trace of its product with the density matrix. Obviously, 
({...) s)b — ((■■■) b)s- The expression (16) is only well defined when all the traces yield final 
results, which should be usually the case for finite-temperature thermal initial bath density 
matrix 9, or for projection-operator 9. More sophisticated short-time approximations have 
been proposed; 6 see below. 

For an instructive example, let us consider next a rather general model of the two-level 
system interacting with a boson- mode reservoir. 12 As usual, to evaluate the density operator 
dynamics, we have to use some approximations. For short time scales, we apply the recently 
developed short-time approximation, 6 rather than the straightforward perturbative expansion 
(15). The Hamiltonian of the system has the form, 

H = H s + H B + V = -ttct z + J2 UJ ka{a k + a x ^2(g k al + g* k a k ). (17) 

1 k k 

Here Hs = —{Vt/2)a z is the Hamiltonian of the quantum system, Hb = J2<^k% a k is the 

k 

Hamiltonian of the bath modes, and the remaining term, V, is the interaction Hamiltonian; 
a k , a\ are the boson annihilation and creation operators; <7j (i = x, y, z) are the Pauli matrices; 
Q > is the energy gap. The eigenstates of a x will be denoted by 

<7 X |±)=±|±), (18) 

where |±) = (|0) ± |1)) /\/2, and |0) =|t), |1) =11) are the ground and excited states respec- 
tively. 

Now, to obtain the short time dynamics of the system, one can use the approximation, 6 



Pi 9 = 0,1 



Pmn{t) = £ {e^ +E "-^- E ™W 2 (m|A*} (//b)(?|^)(^|n)p w (0)e- B2 W^- J? ^ 2 / 4+iC W^-^)}, 

(19) 

^=8£^sin 2 ^coth^, (20) 



C(t) = J2-^r(^kt-smcu k t). (21) 



|2 



k 



Here the Roman-labeled states, are the eigenstates of H$ corresponding to the eigenvalues 
E t = (-l)( i+1 )fi/2, with % = m,n,p,q = 0,1. The Greek-labeled states, \(), are the ± 
eigenstates of a x in (18), with eigenvalues rj^ = (, where ( = /x, v = ±1. Evaluation of (19) 
yields the following expressions, 



Pn(t) = - 



1 + e 



-B 2 {t) 



Pn(0) + - 



-B 2 (t) 



Poo(0), 



(22) 



Pio(*)= 2 e [ 1 + e 



-B 2 (t) 



Pio(0) + - 



3 -i? 2 (i) 



Poi(0). 



(23) 



It transpires that the time-dependence of the density matrix elements within this approxi- 
mation is not exponential. For estimation of the departure from the initial pure state, we 
obtain 

s(t) = 1 - Tr (p 2 ) = \ [l - e~ 2B ^\ {[p n (0) - p 00 (0)] 2 + 4 |p 01 (0)| 2 sin 2 [(fi/2)t - 7o ]} , (24) 



where p i = |poi|e no . 

Note that this result depends on the spectral function B 2 (t), defined in (20). This function 
is obtained by integration over the bath mode frequencies. When the summation in (20) is 
converted to integration in the limit of infinite number of the bath modes, 7 ' 13,14 

B 2 (t) = 8 [ du;N(uj) \g(uj)\ 2 uj- 2 sin 2 ^ coth ^, (25) 



where N(u) is the density of states. In many realistic models of the bath, the density of states 
increases as a power of uo for small frequencies and has a cutoff at large frequencies (Debye 
cutoff in the case of a phonon bath). Therefore, approximately setting 

N(u) \g(tu)\ 2 oc cu n exp (-u/u c ) (26) 



can yield a good qualitative estimate of the relaxation behavior. 13 ' 14 For a popular case of 
Ohmic dissipation, 12 n — 1, the time dependence of B 2 (t) is sketched in Figure 2. One can 
identify the initial stage of quadratic growth of the spectral function, the intermediate region 
of logarithmic growth, and the linear large-time behavior. 

The result (24) illustrates two common features of the decoherence measures. The first 
is the explicit dependence on the initial density matrix elements. The second is the time- 
dependence that involves the spectral frequencies of the system. We note that the spectral 
function B 2 (t) is monotonic. However, the expression (24) also contains oscillatory time 
dependence with the period 27r/Q, which obviously reflects the property of the system's internal 
dynamics rather than its deviation from a pure state. If the initial density matrix is diagonal, 
poi = 0, then we have s(t) = 2 _1 {1 — [poo(0) — pn(0)] 2 e _2B2 ^}, which has no such oscillations. 



Figure 2: Schematic plot of the spectral function B 2 (t) for Ohmic dissipation. 14 ' 7 The "quiet" 
regime t < uj~ x corresponds to B 2 {t) oc (uj c t) 2 , the "quantum" regime u^ 1 <t<(3 corresponds 
to B 2 {t) oc ha(u c t), and "thermal" regime t ^> /3, corresponds to B 2 {t) oc t/(5. 



4 Fidelity 

Writing the total Hamiltonian as usual, 

H = H S + H B + V, (27) 

let us now define the fidelity, 15,16 

F(t) = Tr s [p USaal (t)p(t)], (28) 

where the trace is over the system degrees of freedom, and Pidcai(^) represents the pure-state 
evolution of the system under the constant H$ only, without interaction with the environment 



(V = 0), 



AdeaiW = e- iHat p(0) e l 



(29) 



The fidelity provides a certain measure of decoherence in terms of the difference between 
the "real," environmentally influenced, p(t), evolution and the "free" evolution, pideai(^)- It 
will attain its maximal value, 1, only provided p(t) = pideai(0- This property relies on the fact 
the pideai(£) remains a projection operator (pure state) for all times t > 0. 

Let us consider the two-level system decaying from the excited to ground state, (9,10). In 
this case, there is no internal system dynamics, 

PidealW = ( 1 J , (30) 



^) = ( 1 e -rt). (31) 

and the fidelity is a monotonic function of time, 

F(t) = e' rt . (32) 



Note that the requirement that Pideai(^) is pure-state (projection operator), excludes any 
T > thermalized state as the initial system state. For example, let us consider the infinite- 
temperature initial state of our two level system. We have 

P(0) = PidealW = ( y^, (33) 

which is not a projection operator. The spontaneuos-decay density matrix is then 
The fidelity is constant, 

F(t) = 1/2, (35) 

and it does not provide any useful information of the time-dependence of the decay process. 
We can derive the short-time perturbative expansion 15 for F(t), 

F(t) = l-(t/T F ) 2 + o(t 2 ), (36) 



1/(2t 2 f ) = ((V 2 ) s )b-(((V)s) 2 )b, 



(37) 



which should be compared to (16). 



Let us now turn to the example of the two-level system in the short-time approximation; 
see (17), et seq. Coherent evolution of the two-level system due to the Hamiltonian Hs, is 
described by 

/ Poo(0) poi(0)e <nt \ 



Pideal(^) 



Pio(0)e 



-tilt 



Pn(0) 



(38) 



A straightforward calculation gives 

F(t) = 1 - \ [l - e- B2 ^\ {[ Pll (0) - poo(0)] 2 + 4 |p 01 (0)| 2 sin 2 [(fi/2)t - 7o ]} . (39) 



This result should be compared to (24) While not identical, it has similar features, including 
the system-frequency oscillations. 



5 Norm of deviation 



In this section we propose to use the operator norms 17 that measure the deviation of the 
system from the ideal state, to quantify the degree of decoherence. Such measures do not 
require the initial density matrix to be pure-state. We define the deviation according to 



<j(t) = p(t) - Pideal(*)- 



(40) 



We can use, for instance, the eigenvalue norm, 

||cr|| A = max |Aj| , 



(41) 



or the trace norm, 



iTr 



(42) 



etc., where Aj are the eigenvalues of the deviation operator, (40). 

For more precise definitions, let us consider an arbitrary linear operator A. One of the 
possible ways to define the norm of A is 17 



= sup 



(<p\AA\<p ) 



n l/2 



(43) 



Since density operators are bounded, their norms, as well the norm of the deviation, can be 
always evaluated. Furthermore, since the density operators are Hermitean, this definition 
obviously reduces to the eigenvalue norm (41). We also note that \\A\\ =0 implies that A = 0. 

The calculation of these norms is sometimes simplified by the observation that a(t) is 
traceless. Specifically, for two-level systems, we get 

\W\\x = vVoo| 2 + I coi 1 2 = ^ INI-Tr • (44) 

Therefore, for tow-state systems, we will only consider the eigenvalue norm. For our example 
of the two-level system undergoing spontaneous decay, the norm is 

H A = l-e- rt . (45) 



Thus, in this case ||cr|| A = 1 — F(t). 



Next, consider the two-level system coupled to a bath of modes, in the short-time approx- 
imation; see (17), et seq. The coherent (ideal) evolution of the system due to the Hamiltonian 
H s , is described by (38). One can then obtain the result 



-B 2 {t) 



{[Pn(0) - Poo(0)] 2 + 4 |p 01 (0)| 2 sin 2 [(fi/2)t - 7o ]} 



1/2 



(46) 



which should be compared to (24) and (39). At t — 0, the value of the norm is equal to 0, and 
then it increases to positive values, with superimposed modulation at the system's energy-gap 
frequency. 



6 Arbitrary Initial States, Multiqubit Systems 

The measures considered in the preceding sections quantify decoherence of a system pro- 
vided its initial state is given. However, this is not always the case. Usually, it will be necessary 
to obtain an upper-bound estimate of decoherence for an arbitrary initial state. For exam- 
ple, in quantum teleportation schemes, the quantum state to be teleported is not known. In 
quantum computing, the ideal state of a quantum register is theoretically specified for each 
step of a quantum algorithm. But from the practical point of view, the bookkeeping becomes 
intractable even for few-gate algorithms. Furthermore, even the preparation of the initial state 
can introduce noise. 

To characterize decoherence for an arbitrary initial state, pure or mixed, we propose to 
use the maximal norm, D, which is determined as an operator norm maximized over all initial 



density matrices. For instance, we can define 



D(t) = sup{ |k(t,p(0))|| A 

p(0) 



(47) 



For the two-level system coupled to a bosonic bath, with the dynamics described by the 
Hamiltonian (17), the expression of the maximal norm is elegant and compact, 



D(t) 



1 - e" B2 W 



(48) 



The result is monotonic and contains no oscillations due to the internal system dynamics, as 
shown in Figure 3. 



In principle, the calculation of the maximal norm for a large-scale multiqubit system should 
be formidable task, since it implies maximization over an exponentially large (in the number 
of qubits) set of coefficients of the initial density matrix. However, for specific applications, 
we can use approximate upper bounds calculated by evaluating norms for small subsystems. 
Consider a multiqubit system S consisting of two separate unentangled subsystems S2 and S2, 
at time t, with decoherence norms D\ and D 2 , respectively. Let us denote the density matrix 
of the full system and its deviation as p and a, respectively, and use the same notation with 
indices 1 and 2 for the two subsystems. For brevity, we use the superscript (i) to denote the 
"ideal" density matrices. The overall norm D(t) can then be written as 



D 



sup 

p(0) 



J = sup ( 

p(0) v 



p- P K 



sup 

p(0) 



P1P2 



P1P2 



) = sup ( 

p(0) v 



, (») 



J (49) 



and estimated by using the following sequence of inequalities, 



D 



sup 

p(0) 



&1P2 + Pl 02 



) < SUp (|KP2|| A ) + SUp ( p i \ ) 02 ) , 
K/ P(0) p(0) V A/ 



(50) 



sup 

P(0) 



\WlP2\ 



a) + sup ( 

p(o) v 



(i) 
Pi °2 



) < sup 

x/ Pi(0) 



0^1 



I a) + sup 

P2(0) 



0^2 



(51) 



The inequality D < D\ + D 2 is, of course, approximate, because we have assumed that the 
interaction with the bath, and with each other, left the subsystems unentangled at time t. 
However, it provides an indication that, perhaps with a smart choice of a decomposition 
into subsystems which are initially not significantly entangled, the maximal norm can be 
approximately considered an extensive quantity, linear in the number of qubits, at least for 
short times. It is expected that for larger times, the relaxation process rates are the additive 




Figure 3: Norm of the deviation from the ideal dynamics for the density matrix of a spin 
interacting with an Ohmic bath of bosonic modes, in the short-time approximation. The 
upper curve: the maximal norm D(t); the lower curves: the norms for three different 

initial states p(0), chosen to illustrate the overall pattern. 



quantities growing linearly with the number of qubits, as long as no quantum error correction 
is involved. 

In conclusion, we have considered several approaches to quantify decoherence: relaxation 
times, entropy and fidelity measures, and norms of deviation. The latter measures offer certain 
advantages, and we defined the maximal measure that is not dependent on the initial state 
and, at least for short times, is approximately extensive in the number of qubits. 
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